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The structural and electronic properties of new structural cubic (GaN)i/(ZnO)i super- 
O \ lattice have been investigated using two different theoretical techniques: the full potential- 
linearized augmented plane wave (FP-LAPW) method and the linear combination of local- 
ized pseudo atomic orbital (LCPAO). The new modified Becke- Johnson (mBJ) exchange 
potential is chosen to improve the bandgap of the superlattice and effective masses. The 
bandgap is found to be slightly indirect and reduced from those of pure GaN and ZnO. The 
origin of this reduction is attributed to the p — d repulsion of the Zn-N interface and the 
presence of the O p electron. The electron effective mass is found to be isotropic. Good 
'^ ■ agreement is obtained between two used methods and with available theoretical and exper- 
Q . imental data. 
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.• . 1. Introduction 

'NT 

2^ ! The (GaN)i/(ZnO)i pseudobinary semiconductor alloys, which is discovered in this last 

years [l|, |2| represents a new class of alloys typed as III-V/II-VI. This latter is different from 
classical isovalent semiconductors alloys (i.e, III-V/III-V, II-VI/II-VI). From technological 

'k>( ' aspect, the growth of ZnO on GaN heterostructure was realized using the vapor cooling 
?H ■ condensation system and the related hetero junctions light-emitting diodes LEDs are then 

- - -' fabricated J3|. But the main technological application is the hydrogen generation from water 
photosplitting through the photo-electrochemical cell (PEC). In fact, under visible light, it 
has a much higher water splitting efficiency compared to other oxides |^. The spectac- 
ular feature of this characteristic is attributed partly to the lower bandgap of alloys and 
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in particular of the heterostructure compared to its parent elements GaN and ZnO. Few 
theoretical and experimental studies estimating the optical bowing parameter have been 
reported p, la, l7| . Large difference in the estimation values was than observed. The recent 
calculation of S. Wang and L-W Wang [7| shows that the ordering plays an important role 
in the reduction of the bowing parameter. 

In the other hand, the tendency in this decade of several research groups to elaborate the 
GaN and ZnO in the zinc-blend structure for technological interest, motivate us to contribute 
by a predictive study for this (GaN)i/(ZnO)i superlattice in this structural phase. In fact 
the lattice constants of GaN and ZnO are nearly identical and therefore this superlattice 
in the zinc-blend structure is then possible to growth. Regarding the electronic properties 
from theoretical point of view, the problem of the underestimation of the bandgap using the 
local density and generalized gradient approximations (LDA and GGA) of density functional 
theory (DFT) is very remarkably for ZnO in particular. The recent works have predicted 
the values of 0.710 eV J8| for ZnO and 1.811 eV [9| for GaN which are lowered from the 
experimental ones by 78% and 45% respectively. Consequently, this error affects consider- 
ably the value of the large bandgap discontinuity at the heterointerface for (GaN)i/(ZnO)i 
layered structure. In order to overcome this problem and improve the bandgap of the su- 
perlattice we used the most recent proposed approximation called modified Becke- Johnson 



(mBJ) exchange potential [10|. In fact, this technique is capable to describe with high accu- 
racy the electronic structure of semiconductors and insulators giving rise to a significantly 
improved bandgap values to be much closer to the experimental ones. Because of the lake 
of data about the structural and electronic properties for the studied superlattice we used 
two different ab initio techniques FP-LAPW and LCPAO in order to help understanding 
the related properties and consolidate our results. 

The paper is organized as follows: In section [21 we describe briefly the computational tech- 
niques used in this study. In section [31 we present and discuss our obtained results for the 
electronic properties of the studied superlattice. We summarize our main conclusions in 
section [H 

2. Computational details 

The simulation of the (GaN)i/(ZnO)i structure in the (001) direction is achieved with 
cubic 1x1x1 supercell containing 8 atoms where the atomic positions are those of cubic zinc- 
blende (Figure [T]). The calculations were performed using both full potential-linearized aug- 
mented plane wave (FP-LAPW) and linear combination of localized pseudo atomic orbital 
(LCPAO) methods, based on the density functional theory (DFT) (ill. Il2|. The description 
of the specific parameters to each technique is given as : 

2.1. FP-LAPW 

We used the full-potential linear augmented plane-wave method (FP-LAPW) imple- 



mented in the Wien2k code [13| which self-consistently, finds the eingenvalues and eingen- 



function of the Kohn-Sham equations for the system [ll|, [12|. We used the Generalised Gra- 



dient Approximation (GGA) as parameterized by Perdew, Burke and Ernzenhorf IJ], which 
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Figure 1: The 1x1x1 supcrcell of (GaN)i/(ZnO)i (atomic color map: green for Ga, blue for N, grey for 
Zn and red for O) 



includes the second order gradient components, and Local Density Approximation (LDA) 
both with modified Becke- Johnson exchange potential (mBJ). The core states of Ga, N, Zn 
and O atoms are self-consistently treated and relativistically-relaxed in a spherical approx- 
imation, whereas the valence states are treated self-consistently within the semi-relativistic 
approximation (spinorbit coupling excluded). The valence electron configurations used in 
the calculations are: Ga{3d^Hs^Ap^), N{2s^2p^), Zn{3d^Hs^) and 0{2s^2p^). The wave 
function, charge density and potential are expanded by spherical harmonic functions inside 
non-overlapping spheres surrounding the atomic sites (muffin-tin spheres) and by a plane- 
wave basis set in the remaining space of the unit cell (interstitial region). The maximum 
/ quantum number for the wave function expansion inside atomic spheres is confined to 
Imax = 10. The charge density is Fourier-expanded up to Gmax = 8.5 {RyY^'^. The con- 
vergence parameter Rurkmax which controls the size of the basis set in these calculations 
is set to 7 {kmax is the maximum modulus for the reciprocal lattice vector and Rmt is the 
average radius of the muffin tin spheres). The muffin tin radius values for Ga, N, Zn and O 
were chosen at 1.96, 1.68, 2 and 1.73 atomic units (a.u.) respectively. The reciprocal space 
is sampled by a 12 x 12 x 12 Monkhorst-Pack mesh [l5| with 868/c- vectors in the irreducible 
Brillouin zone. The iteration process is repeated until the calculated total energy of the 
crystal converges to less than 10~^ Ry. 

2.2. LCPAO 

For LCPAO method 
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we used 6x6x6/;; points for the /c-space integration. The 
primitive pseudo atomic basis sets Ga7.0-s3p3(i3/l, Zn6.0-s3p3(i3/l, N6.0-s3p3(i2 and 07. 0- 
s3p3(i3 are used for Ga, Zn, N and O atoms, respectively, where the abbreviation, for exem- 
ple, Ga7.0-s3p3(i3/l, represents the employement of three primitive s orbitals, three primi- 
tive p orbitals, three primitive d orbitals and one primitive / orbitals of a Ga atom which are 
generated with a confinement radius of 7.0 hohr. The standard DFT PP for Zn atom treat 
3(i^°, 4s^ as the valence electrons. In this work, we include the 3s and 3p states as the valence 



electrons to improve accurate calculations [17|. The pseudopotentiels used in the present 



calculations were avaible on the openmx web site 18|. The generalized gradient functional 
developed by Perdew et al [IJ] is used to evaluate the exchange-correlation energy. The real 
space grid techniques are used with the cutoff energy of 450 Ry in numerical integrations and 
solution of the Poisson equation using the fast Fourier transformation (FFT) algorithm. The 
equilibrium positions of ions are reached by the the eigenvector following method where the 
approximate Hessian is updated by the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method 
[l8|. The Hellmann-Feynman forces tolerence is 1 mRy/bohr. All total energies were con- 
verged to 2.0 X 10"^° Ry. 



3. Results and discussion 

3.1. Structural properties 

Figure [2] shows the renormalized total energy of (GaN)i/(ZnO)i superlattice performed 
over a set of different cell volumes. Equilibrium structural properties were obtained by fitting 



the calculated total energies versus volume data to the Murnaghan's equation of states [19 



using both FP-LAPW and LCPAO methods within LDA and GGA approximations. The 
calculated equilibrium lattice parameter a, the bulk modulus B, its first pressure derivative 
B and the total energy Eq are summarized in Tabled] compared to their parent elements GaN 
and ZnO. No experimental data are available for the cubic (GaN)i/(ZnO)i superlattice. The 
only experimental data are given by Maida et al. for the (Gao.87Zno.i3)(No.830o.i6) {a = b = 
3.189i, c = 5.18367i, c/a = 1.625) jioj and the (Gao.93Zno.o7)(No.9oOo.io) {a = b = 3.19A, 
c = 5. 1835^4, c/a = 1.625) [21] in the wurtzite phase at the temperature of 299 K. 



Table 1: Lattice constant a(A), bulk modulus B (GPa), pressure derivative of bulk modulus B' and total 
energy Eo(eV) for (GaN)i/(ZnO)i and their parents. FP-LAPWfLDA) and FP-LAPW(GGA) results are 
obtained by Wien2k [H and LCPAO (GGA) by openmx codes ^. 
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Experiment 




Composition 


FP-LAPW(LDA) 


FP-LAPW{GGA) 


LCPAO (GGA) 


Other calculations 


GaN 


a 


4.4607 


4.5121 


4.5406 


4.52-4.55[22,_26,23l 


4.46[24] 4.50[29] 4.55[30| 4.59[26] 
4.58[27] 




B 


199.5557 


175.8859 


193.3780 


190 [25, 26,^ 


202[24] 176[30| 165.59[26] 




B' 


3.9820 


3.6525 


7.4917 




4.32[24] 3.30[30] 4.24[26] 




Eo 


-54307.19 


-54392.89 


-9471.06 






(GaN)i/(ZnO)i 


a 


4.5859 


4.5866 


4.6332 








B 


147.8437 


145.5577 


128.0100 








B' 


4.1421 


4.3393 


4.5090 








Eo 


-210672.83 


-210632.00 


-17749.60 






ZnO 


a 


4.4938 


4.6204 


4.6497 


4.47[3ll 


4.634[35| 4.6329[3r| 4.637[3^ 
4.53[32] 4.489[33] 4.337[34] 




B 


165.4759 


129.7905 


123.0910 




124[36] 129.19[35] 165.9[331 
139.32[37] 166.649[34] 




B' 


3.9386 


4.0956 


5.2790 




4.33[33| 




Eo 


-50838.65 


-50923.39 


-26028.51 
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Figure 2: Renormalizcd total energy vs lattice parameter of (GaN)i/(ZnO)i superlattice and it parent 



binaries with: (a) FP-LAPW(LDA),Kb3FP-LAPW(GGA) and (c) LCPAO(GGA). Solid lines represent the 



fit by Murnaghan EOS. 



3.2. Electronic properties 

In order to give more accurate results for the superlattice, we have calculated the bandgap 
of GaN and ZnO compounds using both mBJ-LDA and mBJ-GGA. As illustrated in Table 
El our bandgap calculations show that the mBJ-LDA is more improved than the mBJ-GGA 



in accordance with the results of Camargo-Martinez and Baquero [38| for GaN and ZnO 



Figure [3] illustrates the band structure of the (GaN)i/(ZnO)i superlattice obtained by both 



FP-LAPW with mBJ-LDA (Figure [3(a)| ) and LCPAO with GGA (Figure [3(b)| ). By com- 
paring to the binary parent elements GaN and ZnO, the bandgap of the heterostructure is 
lowered in amount considerably (see Table [2]). This finding is an indication of the strong 
positive bowing parameter for the GaNZnO quaternary alloys. S Wang and L-W Wang 
[7| show that ordering plays an important role in increasing the bandgap and reducing the 
bowing parameter for these alloys. As an indication, they obtained a value of 4.8 eV for 
wurtzite structure which is smaller than 11.5 eV [6|. The bandgap of the (GaN)i/(ZnO)i 
superlattice of 1.0419 eV is found to be slightly indirect at R point compared to the direct 
bandgap of 1.1710 eV. As consequence, from technological applications, the superlattice can 
be viewed then as approximately direct semiconductor and therefore to be considered as 
optically active material for photonic and optoelectronic devices. 

The partial densities of states (PDOS) are illustrated in Figure ID The lowest band for 

(GaN)i/(ZnO)i superlattice at 21 eV for mBJ-LDA (wien2k) and at ~ -16 eV for 

LCPAO (GGA) (openmx), is attributed essentially to the O 2s state and the following bands 
of the width 2.857 eV (mBJ-LDA) and 2.966 eV (LCPAO(GGA)) are mainly consisting of 
the localized Ga 3d and N 2s states. The upper valence band (UVB) of the width 7.683 eV 
(mBJ-LDA) and 7.535 eV (LCPAO(GGA)) is formed by the Ga 4s, N 2p, Zn 3d and O 2p 
states, but the maximum of the (UVB) is domination by the N 2p and the Zn 3d states in 
the Zn-N interface of our cubic superlattice. This gives rise to the p — d repulsion which 
has been intensely debated as the responsible of the bandgap narrowing in the wurtzite like 



structure [39|. In the other hand, the non negligible O 2p state in the (UVB) must be also 
considered as the origin of this visible light absorption mechanism produced by this bandgap 
narrowing. This was explained in the basis of the recent photoluminescence (PL) studies by 
some local inhomogeneity of the Zn and O atom densities leading to the impurity levels into 
the bandgap (ie empty impurity just above the valence band and/or filled impurity levels 
just below the conduction band) [40!]. The conduction band originates from the contribution 
of the Ga, Zn, N and O atoms with variable amounts. 

The effective mass of the (GaN)i/(ZnO)i superlattice is an important property to be 
exploited in order to determine the efficiency of its PEC cell. In this context we have 
evaluated the effective masses of electron for the (GaN)i/(ZnO)i superlattice as well as for 
the pure GaN and ZnO, from P point along the principle directions of the Brillioun zone 
using the following expression: 

2m* 



Our results computed by mB J-LDA are given in Table HI showing that the effective mass 
of electron is isotropic for the (GaN)i/(ZnO)i superlattice and equal to 0.198mo. 
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Figure 3: Band structures and total DOS of (GaN)i/(ZnO)i superlattice obtained by: (a) mB J-LDA and 
[(b)]LCPAO(GGA). 



Table 2: Calculated bandgaps by FP-LAPW(niBJ-LDA and mBJ-GGA) with mBJ obtained by Wien2k 
and LCPAO(GGA) by opcnmx. All energies are in eV. 









Present work 


Type 


Experiment 




Composition 


mBJ-LDA 


mBJ-GGA 


LCPAO(GGA) 


Other calculations 


GaN 


Egir" -^ r'^) 


3.1798 


2.9524 


1.6000 


D 


3.25 at RT [41| 
3.30 at 1.6K [4l| 


1.930[24| 1.747[42| 1.811[9] 
1.520[30| 1.706[431 3.13|38| 


(GaN)i/(ZnO)i 


Eg(r" -f r=) 


1.0419 
1.1710 




0.2409 
0.3523 


I 




3.360[40| (wurtzite phase) 


ZnO 


EgCr" ^r'^) 


2.5935 


1.2785 


0.8037 


D 


3.27 at RT[3l| 


0.710[8| 0.641[44| 0.695[45] 
2.549[45| 



Table 3: Calculated upper valence bandwith (UVBW) and total valence bandwith (TVBW). All energies 
are in eV. 





UVBW 


TVBW 


Composition 


mBJ-LDA LCPAO(GGA) 


mBJ-LDA LCPAO(GGA) 


GaN 

(GaN)i/(ZnO)i 
ZnO 


6.7043 6.9739 
7.6828 7.5351 
4.0370 6.0120 


16.6365 15.8717 
21.0845 18.7575 
19.7857 17.7154 



Table 4: The calculated effective mass of electron along each directions from T point obtained by mBJ-LDA. 





Effective mass 


(mVmo) 


Directions 


GaN 


(GaN)i/(Zi 


iO)i 


ZnO 


T ^X 


0.147 


0.199 




0.190 


T^ L 


0.178 


— 




0.197 


T-^ R 


— 


0.181 




— 


T^ M 


— 


0.215 




— 


Average 


0.162 


0.198 




0.194 


Other Gale. 


0.165 [46j 


— 




0.193[47] 
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Figure 4: Partial DOS of (GaN)i/(ZnO)i superlattice for : (a) Ga, (b) N, (c) Zn and (d) O calculated with 
mBJ-LDA. 
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Figure 5: Partial DOS of (GaN)i/(ZnO)i superlattice for : (a) Ga, (b) N, (c) Zn and (d) O calculated with 
LCPAO(GGA). 



4. Conclusion 

In this study we have investigated the structural and electronic properties of new struc- 
tural cubic (GaN)i/(ZnO)i superlattice. The structural optimisation as well as the electronic 
band structure are obtained using both full potential-linearized augmented plane wave (FP- 
LAPW) and linear combination of localized pseudo atomic orbital (LCPAO) methods. The 
bandgap of the parents elements GaN and ZnO are obtained to be more improved with 
mBJ-LDA than with mBJ-GGA. Our estimated mBJ-LDA bandgap value of the studied 
superlattice is slightly indirect (1.042 eV). This later is lowered considerably from those of 
the binary parent elements GaN and ZnO, indicating to the strong bowing parameter as is 
it observed in previous studies for the wurtzite structure. We have found that the origin of 
this decreasing is attributed to the p — d repulsion of the Zn-N interface and the presence 
of the O p electron. The electron effective mass is found to be isotropic. Good agreement is 
obtained between two used methods and with available theoretical and experimental data. 
We expected that our predictive results of the cubic (GaN)i/(ZnO)i superlattice can be 
used as references for futur theoretical and experimental works. 
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